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1. Introduction 

Some of the most stringent tests of the standard model and most promising windows into 
possible new physics come from process which occur at second order in the Fermi constant Gp. 
These include the standard model prediction for indirect CP violation in the K° — Z° system and 
the bound on interactions which change strangeness by two units coming from the small Kp — K$ 
mass difference. Such 0(G F ) amplitudes involve the exchange of two intermediate bosons which 
requires that each pair of vertices joined by such a W ± exchange are separated by a short distance 
on the order of 1 jmyj. However, the spatial separation between these two pairs of vertices need not 
be small, and "long distance" separations of the order of 1/Aqcd are possible. 

This is well illustrated by the familiar Wigner-Weisskopf formula describing the time evolution 
of the K°-K system: 




M 00 M oT) \ _ i_ ( r 00 r o0 
M w M m) 2 V r ooroo 




(1.1) 



where the 2x2 matrices M and Y are given by: 

r i7 = 2rcW dE(i\H w \a(E)){a(E)\H w \j)8(E-m K ) (1.2) 
M, = y> f dE Ww\a(E))(a(E)\Hw\j) , (L3) 

We are using the subscripts and to represent the K° and Z° states. The generalized sum over 
a and integral over the energy represents the sum over a complete set of energy eigenstates nor- 
malized as (a'(E')\a(E)) = 8{E' — E)8 a / a . These two matrices determine the two important 
quantities, mg s — mx L and Ek'- 

m Ks -m KL = 2Re{M (,} (1.4) 
i_ flmM 05 -ilmT 05 ) IrnAp 
2 |ReM o-^Rer oJ ReA 

The absorptive part, given by the energy conserving matrix r i; , is formed from products of 
first-order weak amplitudes which are now the targets of increasingly precise lattice calculations, 
see for example Refs. [jl|] and The dispersive part, M, 7 - is intrinsically second order in Gp and 
contains both a short-distance part were both exchanged IV's are separated by distance scales much 
smaller than 1 /Aqcd and the long-distance part discussed above. 

The computation of the short distance contributions to Sk or m# s — mk l is now a highly suc- 
cessful application of lattice QCD. At the scale of QCD, phenomena taking place at a distance scale 
of 1 /mw can be accurately represented by an effective four-Fermi operator whose coefficient can be 
computed in electro-weak perturbation theory and whose matrix element between kaon states can 
be computed using lattice QCD. For Ek these short-distances dominate and the long-distance part is 
estimated to be a few percent correction []3p. However, current results for this short contribution to 
Ek are now accurate on the 5% level, giving the unknown long-distance part increasing importance. 
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For the Kl — K$ mass difference the size of the long-distance contribution is less certain and could 
provide as much as 50% of the full result. 

It is the lattice calculation of these long distance parts which is the subject of this talk. We will 
briefly review the approach presented in Ref. [Q] for the calculation of the long-distance contribu- 
tion to the Kl — K$ mass difference and then describe its generalization to the more complex case 
of e K - 



2. Strategy for lattice calculation Kl — K$ mass difference. 

As presented in Ref. [Q] a Euclidean-space, finite volume calculation of the Kl — K$ mass dif- 
ference Amjc contains three important ingredients. The first is a relation between the second-order, 



infinite- volume, principal part integral in Eq. [L3] giving Amx = 2Re(Mg ) and a finite- volume, dis- 



crete perturbation theory sum. This relation is obtained through a generalization of the first-order 
formula of Lellouch and Luscher for the K — > kk decay width [§J. The second ingredient is a 
computational strategy for evaluating that perturbation theory sum using lattice QCD. While with 
sufficient statistics the resulting lattice calculation will correctly determine the long-distance part 
of AniK, the short distance contribution to this result will be incorrect, reflecting the details of the 
lattice regulator. Thus, the third ingredient is a method to replace this erroneous short-distance con- 
tribution with the correct short distance part of the continuum theory. We will now briefly review 
each of these in turn. 

The relation between Am^ and a finite volume perturbation theory sum follows from the 
Luscher finite-volume quantization condition [||] relating the allowed finite volume energies, E n , 
for a two-pion system and the % — % scatting phase shift: 

Q(k n L/2K) + 8o(E n ) + 8 w (E n )=nn (2.1) 



where k n = ^jE} x jA — m\ and the known function (j)(y) is defined in Ref. [g]. Here we have divided 
the k — n scatting phase shift into the strong s-wave phase shift 5q and that caused by the resonant 
scattering through the weakly coupled K$ state. Since the location of the kaon pole in 5w(E n ) is 



shifted by Amk, Eq. [2JJ can be used to relate Amk to the finite-volume, second-order perturba- 
tion theory sum which determines E n . It is this argument which we will generalize to determine 
Im(M 5() ). 

Since the resulting finite volume sum that must be evaluated to determine AniK is closely 
related to perturbation theory, the entire sum can be obtained from an integrated Green's function 
constructed from four operators: two interpolating operators which create an initial K° state at 
a time ti and destroy a final K state at tf and two AS = 1 weak operators Hy/{t\) and Hwfo) 
evaluated at intermediate times and integrated over a range tf^>tj,> 4 > t a S> ?, for k = 1, 2. The 
term proportional to length of the integration interval, % — t a gives the desired perturbation theory 
sum which determines Amk- While the evaluation of this time-integrated four-point function poses 
substantial computational challenges, our first exploratory study gives encouraging results 

Since ArtiK receives both short- and long-distance contributions, the integrated amplitude de- 
scribed in the previous paragraph will receive important contributions from the region in which 
the two weak vertices coincide. If a charm quark is included to exploit GIM cancelation a sub- 
stantial, incorrect short-distance contribution will remain, of size hv{mwa). Fortunately, the same 
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techniques which allow this continuum short-distance part to be precisely defined and evaluated in 
a lattice calculation can be employed here to isolate and replace this lattice artifact by the correct 
continuum piece. If the four quarks making up the two kaon interpolating operators in the above 
four-point function are instead treated separately and the resulting Green's function evaluated in 
Landau gauge and given large, off-shell momentum the short distance contribution can be evalu- 
ated numerically in a fashion that is free of infrared singularity allowing a simple substitution by 
the correct continuum short distance part. Again, first numerical experiments [f7|] suggest that this 
can be done without difficulty. 



3. Extension to £k 

We now turn to the central topic of this talk, the extension of the method reviewed above 
to the case of Im(MQ ) and Ek- Since there are no apparent added difficulties associated with the 
numerical evaluation of the resulting finite- volume expression or the correction of its short-distance 
part these topics will not be discussed further. However, on first-sight using a Lellouch-Luscher 
style argument to obtain a finite volume formula for Im(M 00 -) may seem unlikely to succeed. First 
e K describes the mixing of the K s and K L states and the K L state with its predominate decay into 
three pions would appear inaccessible to Luscher's formalism which applies only to two-particle 
scattering. Further the energies that are the basis of the Lellouch-Luscher method are intrinsically 
CP conserving suggesting no sensitivity to the CP violating phase of M q. 

Fortunately both of these problems can be avoided by adding a fictitious, AS = 2 "superweak" 
term to the Hamiltonian whose effects can be combined with (but will not alter at leading order) 
those of the standard weak interactions: 

H sw = (cOr + imi) {s(\ + Y') r fd 1(1 + 7 5 )y /i <i} +hermitian conjugate. (3.1) 

The quantities co r and ft), are then chosen so that 

ev(*° 1 LL \K°) + ReM 5o = (3.2) 
^•(iVttlZ ) » |r»|, (3.3) 



for all i and j. Here Oll is the operator within the curly brackets in Eq. [3J. With the off-diagonal, 
real part of Mjj canceled and the off-diagonal term proportional to ft), dominating the off-diagonal 
parts of r ;y , the eigenstates of the time development operator in Eq. [Tj] become the unfamiliar 
combinations K± = (| K°) ± i\K ) )/ V2 — both states which decay predominantly into two pions, 
consistent with Luscher's finite volume formalism. In addition, the large contribution of ft), insures 
that the two % — % resonances associated with the states K± are non-overlapping and will each 
contribute independent resonant behavior to that scattering. The eigenvalues of the 2x2 time 



development matrix on the right-hand side of Eq. |Ll| for the states \K±) can be written: 



l - (M m +M m ) T ^{K^OllIK 1 ) +ImM o 



+4( r °° + r oo^' r oo±<r o) (34) 



Moo T ((Oi(K Q \ LL \K°) + ImM o5 ) + i~ (r 00 =F Imr oo -) (3.5) 



4 



Long-distance contribution to £k 



Norman Christ 



We first examine the finite volume problem. Here the box size L is adjusted so that we have 
three nearly degenerate states: the % — % state |ran) with energy E no as well as the two single-particle 
states, \K°) and \K°). Using second order, degenerate perturbation theory, we can determine the 
finite volume energies of interest as the eigenvalues of the 3x3 matrix: 



m K + 



E \{n\H w \K°)\- 
n=in m K -E„ 



E (T\H w \n)(n\H w \K°) | g m „ 



E {K°\H w \n)(n\H w \K ) g 
n^no rn K -E„ _ ~ r 

E l^l^lg )! 2 



V 



m K -E„ 

(n \H w \K°) 



m K + 



n^no m K -E„ 
(n \H w \K°) 



(K \H w \n ) 

E \(n\H w \n )\ 2 
n^K ^ E„ Q -E„ J 



(3.6) 



with the rows and columns corresponding to \K°), \K") and |no) in increasing order. Here the 
complex quantity dm is given by: 



5m = (co r + /w,-)(^ |C?LL|^' 



(3.7) 



Using the CP and time-reversal operations, the 3x3 matrix in Eq. 3.6 can be written as: 

^ E d \ ^ e 

Jt* E di e- iS °£. 
^e iS °£/* e iS °£/ E d2 



(3.8) 



The eigenvalues of this matrix can be determined by diagonalizing the large part of the matrix 
proportional to the first-order amplitude srf and then applying perturbation theory to diagonalize 



the remainder to order 



Eq 



Ed\ 
1 



Re{^(=c/*) 2 } 



{E dl +E d2 )±V2\#/\ + 



Re 



* \21 



2\*/\< 



(3.9) 
(3.10) 



We now require that the Luscher finite volume condition, Eq. 2.1, be obeyed when the total 
7T — 7r phase shift, written as a combination of the usual strong s-wave phase shift, the resonant 
scattering through the \K±) states and a standard second-order-weak Born term is evaluated at the 
finite-volume energy eigenvalues E± given in Eq. [3.10t 



nn 







' » + So(E±) + arctan + V ± >' + arctan ' 

\ M + — h± J 



2% 



-* I 



M_-£± 
\{P\H w \tck(E))\ 2 

E-Ea 



(3.11) 



Here M± and F± can be determined from the real and imaginary parts of the infinite- volume eigen- 



values given in Eq. ^5 



r± 



m K + Moo T («i 1 LL + ImM o5 

r 00 =Fimr o5 . 



(3.12) 
(3.13) 



If the finite volume is adjusted so that the zero-order % — % energy E no = m^, then to zeroth 
order in Hy/, Eq. 3.11 reduces to the original Luscher quantization condition relating the strong 
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s-wave phase shift 5q(E±) and the finite- volume 71 — 71 energy E„ . If, for the same value of E, 



-<«o> 



Eq. p.ll| is expanded to first order in Hw then the Lellouch-Luscher relation between the finite 
volume K — 7171 matrix element and the A* — >• %% width Too is obtained. 

Considerably more algebra is needed to obtain a relation between the finite volume sum in 



and the infinite volume amplitude M 00 -. As in the earlier study of Anix, we first evaluate Eq. p.l 1 
for a volume which makes the 71 — 71 energy E„ close to on the scale of QCD but sufficiently 
different that a perturbative expansion in Hw/(E„ —m£) is possible. Equation. 3.1 1 can then be 
evaluated in a power series in H^. The resulting formula contains a pole at E llQ = rtiK- Requiring 
that this pole term vanish gives the standard Lellouch-Luscher relation. Requiring that the non-pole 
piece vanishes at E nQ = mg then relates the infinite- volume, second-order Born term in Eq. 3.11 



and the finite- volume, second-order sum in the lower left corner of the matrix in Eq. 3.6. 

In the next step, Eq. 3.1 1 is evaluated for a volume making the two-pion energy E no = and 
the result expanded to second order in Hw, including the resonant denominators. The expression 

■4 
W 



that results contains denominators of order Hy, and numerators of order H% and can be written: 







d(<!> + 8o) 




■Mi 



oo 



"F70v 



(K°\H w \n)(n\H w \K 



m K -E n 



M, 



^±^|< n „|»w|K»>l 




(3.14) 



giving a finite volume expression for a combination of the infinite volume amplitudes Moo and M 00 -. 

To determine M 00 - alone, a second equation is needed. As in the simpler case examined Ref. [EJ], 
a second relation follows is we require that the expectation value of the matrix Mjj for the state 
\K) = (^\K°)-£/*\TC))/(2\^\ 2 ) agree with the corresponding finite volume expression up to 
terms exponentially small in the size of that volume. Here the state \K) is chosen so that it does not 
couple to the two-pion, on-shell state, justifying the neglect of finite volume effects. The resulting 
equation can be combined with Eq. 3.14 to eliminate Moo giving the relation: 



= Re< 



+ 



/*\2 



'°\H w \n)(n\H w \K°) 



m K -E n 



■M, 



(3.15) 



1 



3e 



2 d& * 



<j — \K \H w \no){n Q \M w \K 



dE no \ 



Since the K — > %% amplitude is non-zero with a phase that is independent of the phase of the 
amplitude in curved brackets, (. . .), we can conclude that (...) = 0, our desired finite-volume ex- 
pression for M q. The real part of this equation gives the earlier result for Re{M 00 -}. The imaginary 
part gives the new result for Im{M (,}: 



Im{M o5 } 




Im{(K°\H w \n)(n\H w \K°)} 



m K - E n 



+ 



1 



d(<j>+So) 
5e 



ld 2 (0 + *) 



dE 2 



d f d(<l> + So) » o, 
<j — (A \Hw\no){n \HwK 



dE n \ 



{K Q \H w \n ){n Q \H w K U 



(3.16) 
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Thus, by exploiting the interference between the complex, CP-violating amplitude 80 and the 
second-order, Wigner-Weisskopf mass matrix, we have been able to determine the CP-violating 
part of M q from a finite-volume energy which is intrinsically CP-even. 

Here the left hand side of Eq. 3.16| is the infinite volume quantity which determines the long- 



and short-distance dispersive parts of Ek- The right-hand side contains only finite volume matrix 
elements which can be evaluated in a lattice QCD calculation. The similarity of the results for 
Re{M 00 -}, Im{M () Q} and Moo (not displayed) suggests that a more direct derivation may be possible. 

4. Conclusion 

Given the control over the singularity associated with the principal part and related finite vol- 
ume effects implied by these results and those in Ref. [Q], both the Kl — K$ mass difference Am# 
and the complete, dispersive part of Ek can in principle be computed using lattice methods. In 
fact, the first exploratory results of Jianglei Yu, presented in Ref. 0, suggest that such a lattice 
calculation may be possible with the next generation of high performance computers. 

However, many obstacles must be overcome: (a) This initial work suggests that the short- 
distance parts of AniK and Ek can be properly treated by a single subtraction, defined using Rome- 
Southampton techniques. Such a single subtraction is adequate only if GIM cancelation has been 
realized in the lattice calculation, requiring the inclusion of the charm quark which in turn requires 
a small lattice spacing, (b) Many more operators and contractions must be included in a complete 
calculation than have been attempted in Ref. [0]. (c) These initial calculations use a relatively 
massive, 410 MeV pion. The successful subtraction of an exponentially growing single-pion con- 
tribution demonstrated in Ref. |Q] will become more difficult as the pion mass is reduced. Even 
greater difficulty will result from the vacuum subtraction needed when disconnected graphs are 
included. We believe that these difficulties can be overcome with improved numerical methods 
and the sustained, 100 Tflops capability of the next generation machines now being installed. The 
author acknowledges the important contributions of his RBC/UKQCD collaborators to this work. 
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